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Electron capture decay of 
indium-Ill human carbonic anhydrase I: 
A time differential K X ray coincidence 
perturbed angular correlation study 

Christopher 
Research Computing Facility 
Mayo Foundation, Rochester, MN, USA 

The relaxation effects in the perturbed angular correlation spectra of 111 In human carbonic an- 
hydrase I (HCA I) are the result of chemical transmutation and/or the complex Auger cascades 
that follow the electron capture decay of In. Time differential K X ray coincidence perturbed 
angular correlation (PAC) spectroscopy shows that these relaxation effects are independent of the 
Auger cascade intensity. This suggests that chemical transmutation is responsible for the relaxation 
effects, and that bond breaking and damage product formation around the decay site resulting from 
localized energy deposition by Auger and Coster-Kronig electrons probably occur in the microsec- 
ond time regime. Numerical simulations of chemical transmutation relaxation effects in the time 
differential PAC spectrum of 111 In HCA I are also presented. 



Introduction 



Perturbed angular correlation spectroscopy (PAC) is 
a potential technique for investigating hot atom chem- 
istry following electron capture decay of a radionuclide 
probe 0,0- As is the case with Mossbauer emission 
spectroscopy, it may be possible to identify the chemi- 
cal form of the probe nuclide coordination complex dur- 
ing the nuclear lifetime 0, Q . The average Auger and 
Coster-Kronig electron yield is significantly diminished 
when the vacancy cascade following electron capture in- 
cludes the emission of a K X ray as compared to when 
it does not. It is in principle possible to study the de- 
pendence of hot atom chemistry on the intensity of Auger 
and Coster-Kronig emissions by detecting the nuclear ra- 
diations of perturbed angular correlation or Mossbauer 
emission spectroscopy in coincidence with K shell X 
rays 0, Q . Given the very limited knowledge of hot atom 
chemistry of electron capture nuclides that are possible 
perturbed angular correlation probes, one may ask only 
a yes or no question. Does or does not the hot atom 
chemistry depend electron emission intensity? To answer 
this, a chemical transmutation model and a dose depen- 
dent damage model are constructed for the perturbed 
angular correlation spectra following electron capture. 
Both models encompass several alternative pictures of 
the hot atom chemistry. The simplest hot atom chem- 
istry of the chemical transmutation model is very sim- 
ilar to j3 + decay on the nanosecond time scale. Note 
that j3 + decay is energetically forbidden for m In. The 
multiple Auger ionization processes are followed by com- 
plete electron recombination that leaves the daughter nu- 
clide in a metastable state with both the atomic number 
and valence decreased by one unit from the parent. The 
metastable state is the result of the change in charge 
rather that charge neutralization or Auger electron irra- 
diation. The metastable state of the daughter nuclide 



decays within tens of nanoseconds to the ground state 
coordination complex of the daughter and ligands. How- 
ever, any hot atom chemistry involving an excited state of 
the daughter that decays thereafter into the ground state 
complex is equally consistent with the chemical transmu- 
tation model. The only requirement is that the excited 
state be independent of Auger emission intensity. The 
dose dependent damage model pictures radiolytic frag- 
mentation of probe coordination complex ligands. The 
amplitudes and rates of ligand fluctuations increase with 
the intensity of the Auger emissions. Though uncer- 
tain about the details of hot atom chemistry both mod- 
els predict perturbed angular correlation spectra. The 
comparison of these spectra with experimental data in- 
dicates whether or not the hot atom chemistry during 
the nanoseconds following electron capture depends on 
Auger emission intensity. 

A partial decay scheme for ln In and lllm Cd is shown 
in Fig. ^ Perturbed angular correlation spectra of the 
isomeric m " l Cd decay provide the quadrupole interac- 
tion parameters of the ground state coordination com- 
plex in the chemical transmutation model and the ligand 
relaxation parameters at zero dose in the dose depen- 
dent damage model. Following m In electron capture 
decay the m Cd daughter is in a 120 picosecond nuclear 
excited state that initiates the 171-245 keV gamma cas- 
cade. Processes occurring faster than the lifetime of this 
initial nuclear state, including dissipation of thermal hot 
spots |3j 13 , do not influence perturbed angular correla- 
tion spectra. The spectra reflect the hot atom chemistry 
occurring during the lifetime of the intermediate state. 
Indirect effects involving the diffusion of water radicals 
are not observed. A great variety of radiochemical and 
radiobiological experiments show that direct and indirect 
processes together result in damage that is proportional 
to the Auger emission intensity. For example, when zinc 
bovine carbonic anhydrase is X ray irradiated just above 
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FIG. 1: Partial decay scheme for 111 In and 11 lm Cd 0. The 
half-life of each level is labeled. The 171-245 keV gamma 
cascade follows the electron capture decay of nl In; the 48.5 
minute lllm Cd state decays by the 151-245 keV gamma cas- 
cade. Both cascades share the same 245 keV 5/2+ 111 Cd 
intermediate nuclear state. 



the zinc K edge, enzyme activity and zinc release assays 
indicate significant Auger inactivation 9]. The more in- 
tense Auger emissions accompanying m In decay can be 
expected to inactivate HCA I even more effectively. The 
combined direct and indirect effects are Auger emission 
intensity dependent. An Auger dependent effect on the 
nanosecond time scale can be detected with K X ray co- 
incidence perturbed angular correlation spectroscopy. 



Materials and Methods 

Human carbonic anhydrase I was isolated from 
freshly outdated erythrocytes The apoenzyme 

was obtained by dialysis at 4°C against pyridine-2,6- 
dicarboxylic acid [13 . HCA I activity was assayed spec- 
trophotometrically with p-nitrophenyl acetate (Sigma 
Chemical Co., St. Louis, Missouri) as substrate |12| . 
Conventional time differential gamma-gamma perturbed 
angular correlation spectra of 111 In labeled HCA I were 
measured on a movable two detector perturbed angular 
correlation spectrometer jj^- A 2 by 2 inch cylindrical 
BC-404 plastic scintillator with 0.25 inch centered side 
bore for internal sample mounting and RCA 8575 photo- 
tube detected the K X rays. Coincidence gating electron- 
ics, including a time to pulse height converter for coinci- 
dence timing K X rays and the gamma-gamma cascade, 
selected either K X ray coincidence or anti-coincidence 
spectra for storage and analysis. Carrier free m In in 0.05 
molar HC1 (Du Pont de Nemours & Co. Inc., Wilming- 
ton, Delaware) was ordered at two week intervals. Apo- 



HCA I was incubated for 48 h with m In in 100 mM 
HEPES buffer (Research Organics Inc., Cleveland, Ohio) 
at pH 7.7 with sodium citrate acting as carrier for in- 
dium. In order to hold constant both the volume and 
total indium concentration, the specific activity of the 
m In stock was adjusted with cold indium chloride. Af- 
ter labeling with indium, the effects of HCA I rotational 
motions on the perturbed angular correlation spectrum 
were diminished by adding sucrose to a final concentra- 
tion by weight of 50%. Seven labeled protein samples 
were prepared from each shipment of In. A perturbed 
angular correlation spectrum was accumulated for 46 h 
from each sample. The first and last samples were run 
as incubation controls. These samples were measured for 
an initial twelve h as m In-citrate; HCA I was added, 
and the measurement was continued for another 34 h. 
The time integral perturbation factors accumulated for 
20 min intervals monitored the HCA I labeling reaction. 
All spectra from shipments that did not give excellent la- 
beling were discarded. Three of the five non-control sam- 
ples per shipment were measured in K X ray coincidence 
mode and two in K X ray anti-coincidence mode. A total 
of 52 coincidence mode samples and 36 anti-coincidence 
mode samples were combined for analysis. All sample 
preparations and measurements were done at 21 ± 1°C. 

Electron capture gamma-gamma perturbed angular 
correlation spectra measure the ensemble average elec- 
tric quadrupole relaxation of the daughter nucleus spin. 
Magnetic dipolc interactions are assumed to be negligi- 
ble. The measurable parameters are model dependent. 
All models must contain one or more rate constants that 
characterize the time dependence of the quadrupole in- 
teractions. In the chemical transmutation model, spin 
relaxation is the result of the decay of a metastable 
state of the daughter nucleus coordination complex. This 
metastable state might be a valence electron, coordi- 
nation geometry, or ligand conformation state. In any 
case, the quadrupole interaction parameters of the ini- 
tial metastable state are the frequency mean, frequency 
standard deviation, asymmetry, and the mean lifetime of 
the state. The initial state decays into the ground state 
coordination complex. The interaction parameters of the 
ground state are the frequency mean, frequency standard 
deviation, asymmetry, and Euler rotation angles relative 
to the initial state. The initial and ground state interac- 
tion frequencies are independent Gaussian distributions. 
Each individual daughter nucleus interacts with the ini- 
tial and ground state coordination complexes at fixed 
quadrupole frequencies that are randomly and indepen- 
dently sampled from the respective Gaussian distribu- 
tions. The lifetimes of the initial interactions are given 
by an exponential distribution. In addition to relaxation 
from coordination complex decay, the chemical transmu- 
tation model includes relaxation from overall molecular 
tumbling by multiplying the spectrum with an exponen- 
tial factor. 
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The dose dependent damage model ascribes the spin 
relaxation to ligand reorientation in the coordination 
complex. The quadrupole interaction parameters of the 
coordination complex are the mean frequency, frequency 
standard deviation, asymmetry, rotation angle standard 
deviation, and jump rate. These five parameters depend 
linearly on dose to the coordination complex and have 
a baseline value for the zero dose coordination complex. 
The electric field gradient orientation distribution is de- 
fined by rotations from a reference orientation, where the 
rotation angle distribution is Gaussian and the rotation 
axis distribution is isotropic. Note that though there is a 
single reference orientation in this interaction model, the 
perturbed angular correlation spectrum is for an isotropic 
source because the interactions of each daughter nucleus 
are averaged over random orientations when the spec- 
trum is computed. The interaction frequency distribu- 
tion is a Gaussian. Both the interaction frequency and 
orientation are sampled after every jump. The time in- 
tervals between jumps are exponentially distributed. 

In both the chemical transmutation and dose depen- 
dent damage models, the interaction of each individual 
daughter nucleus is piecewise static and each set of model 
parameters specifies a probability distribution function 
for piecewise static interactions. The perturbation factor 
(which is loosely refered to as the perturbed angular cor- 
relation spectrum), is the average over this probability 
distribution of the perturbation factor of each piecewise 
static interaction. For a piecewise static interaction aver- 
aged over random source orientations, the perturbation 
factors are [Tiij . 

Gkk(t) = {(-l) 2/+m ° +m6 

m a miyN 

(I i k W i i k \ 

\ m' a -m a N J \m' b -m b N J 
x (m b \A(t)\m a )(m' b \A(t)\m' a )*}, (1) 

where k is an even index, t is the delay time between 
emission of the two gamma rays from the daughter nu- 
cleus gamma-gamma cascade, / is the spin of the inter- 
mediate nuclear state, m , m' a , mi,, and m' b index mag- 
netic substates of the intermediate state, the matrix ele- 
ment coefficients are the Wigner 3-j symbols, and A(t) is 
the piecewise static time evolution operator. Since only 
G22 (t) is usually measurable, it is refered to as the pertur- 
bation factor. The index N takes all integer values with 
absolute value less than or equal to k and the primed 
magnetic substate indices are given by m' a — m a + N = 
and m' b — nib + N = 0. Since the perturbation factor 
is always computed at a series of time points, it is ap- 
propriate to express the piecewise static time evolution 
operator as a series operator product of interval evolu- 
tion operators. The interval evolution operator for the 
time interval (t',t) is a series operator product of static 



evolution operators, 

J 

A(t',t) = JJ e - iH ^ /n , (2) 
3= J' 

where Hj is the j th time independent interaction Hamil- 
tonian, J is the index of the first Hamiltonian with ac- 
tive interval overlapping (i', t), J' is the index of the last 
Hamiltonian with active interval overlapping (t',t), and 
tj is the time overlap of (t',t) and the active time in- 
terval of the j th Hamiltonian. Matthias et al. give 
the electric quadrupole Hamiltonian matrix elements for 
arbitrary intermediate spin, interaction frequency, asym- 
metry parameter, and y convention [l5| orientation. The 
quadrupole interaction frequency is specified by the an- 
gular frequency uiq. By definition, loq equals 3 times the 
quadrupole interaction frequency for integer spin and 6 
times the quadrupole interaction frequency for half inte- 
ger spin. When the electric quadrupole Hamiltonian is 
axially symmetric uj^Ti equals the smallest nonvanishing 
eigenvalue difference. The angular frequency ujq is fre- 
quently reported with correct numerical value but with 
the erroneous dimension of Hertz [l6| . 

The piecewise static perturbation factors were evalu- 
ated by repeated extension of the time evolution oper- 
ators in increments of 0.025(27r/u;o), where u>q was the 
mean angular frequency of the ground state or zero dose 
coordination complex. The evolution operator at each 
time point was extended to the next time point by ma- 
trix multiplication with a time interval evolution opera- 
tor. Each time interval evolution operator was a product 
of static evolution operators with the Hamiltonian in- 
dices and overlap intervals indicated in Eq. J2J). The in- 
teraction Hamiltonians were diagonalized |l7j : the static 
evolution operators were expressed in diagonal form with 
the eigenvalues, transformed back into the magnetic sub- 
state representation with the eigenvectors, and multiplied 
together to give the time interval evolution operators. 
Since in the parameter domain of interest the pertur- 
bation factor time intervals were an order of magnitude 
smaller than the typical static Hamiltonian active time 
interval, most interval evolution operators were a single 
static evolution operator and each Hamiltonian diagonal- 
ization generated many interval evolution operators. The 
time interval evolution operators were repeatedly multi- 
plied to give the evolution operators at a series of 160 
time points. The matrix elements of the evolution oper- 
ators were inserted in Eq. Q to give the piecewise static 
perturbation factors at these time points. 

The piecewise static perturbation factor probability 
distribution function was averaged by Monte Carlo in- 
tegration 0]. Perturbation factor error estimates were 
calculated from the deviation of subaverages around the 
overall average as a function of the spectrum delay time. 
Gaussian interaction frequency and rotation angle dis- 
tributions were generated by summing uniform random 
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variates. Excited state and orientation lifetimes were 
generated by logarithmic transform of uniform random 
numbers on the unit interval. Isotropic rotation direc- 
tions were generated by an acceptance-rejection tech- 
nique |l9j . Random numbers uniform on the unit interval 
were generated by a 32-bit linear congrucntial pseudo- 
random generator with divisor 2 31 — 1 and multiplier 7 5 . 

Initially the qualitative dependence of the perturbation 
factor on parameter space of the chemical transmutation 
and the dose dependent damage models was coarsely sur- 
veyed. The chemical transmutation model metastable 
state lifetime, ground state rotation angle, and relative 
magnitude of the metastable and ground state mean in- 
teraction frequency parameters were varied. Based on 
the preliminary survey and known ground state coordi- 
nation complex parameters, a crude estimate was made 
of the parameter point that fit the data. This fit was re- 
fined by varying the unknown parameters one at a time 
until the residuals were about the same magnitude as the 
experimental error. The chemical transmutation model 
ground state frequency mean, frequency standard devia- 
tion, and asymmetry and the overall tumbling relaxation 
time constant were fixed at the known values. The Euler 
rotation angles relative to the initial state were fixed at 
values found in the preliminary survey to give strong spin 
relaxation. The initial metastable state was assumed to 
be axially symmetric, and the frequency mean and fre- 
quency standard deviation were varied. 

The preliminary dose dependent damage model param- 
eter survey fixed the mean frequency, frequency standard 
deviation, and asymmetry parameters at the known zero 
dose coordination complex values, and varied the rota- 
tion angle standard deviation and jump rate parameters. 
By comparing of the results of this search and the per- 
turbation factor measured by Bauer et al. [2CJ with the 
assumption that spin relaxation is due to molecular tum- 
bling, values were identified for the rotation angle stan- 
dard deviation and jump rate parameters approximately 
equivalent to the molecular tumbling time constant of the 
zero dose coordination complex. The preliminary search 
also yielded the initial parameter point for fitting the 
dose dependent damage model to the data. The mean 
frequency, frequency standard deviation, asymmetry, ro- 
tation angle standard deviation, and jump rate parame- 
ters were all varied. 

During the preliminary parameter space searches and 
parameter refinement, the perturbation factors were 
computed by averaging 10 3 piecewise static interactions. 
The standard deviations of these perturbation factors 
were fairly independent of spectrum delay time and typi- 
cally about 0.01. The displayed perturbation factors were 
averaged over 10 4 interactions to give a standard devia- 
tion of about 0.003. All calculations were done on a DEC 
VAX 3600 or Silicon Graphics 4D workstation. 
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20 40 60 80 100 120 
t (ns) 

FIG. 2: Chemical transmutation model for HCA I perturbed 
angular correlation spectra. The triangles with standard er- 
ror markers are the K X ray anti-coincidence spectrum and 
the stars with standard error markers are the K X ray coinci- 
dence spectrum. The solid curve is an approximate fit of the 
chemical transmutation model to the data. The dotted curve 
is the least squares fit of Bauer et al. |H to the lllm Cd HCA 
I perturbed angular correlation spectrum. 



Results and Discussion 

The measured K X ray coincidence and anti- 
coincidence perturbed angular correlation spectra for 
nl In HCA I in 50% sucrose are shown in both Figs. [3 
& |3J Figure [21 shows theoretical spectra for the chem- 
ical transmutation model at parameter values listed in 
Table [I] and Fig. [21 shows theoretical spectra for the dose 
dependent damage model at parameter values listed in 
Table [H] The procedure for approximately fitting the 
models is given in the methods section. Relaxation in 
the chemical transmutation model spectra results from 
metastable coordination complex decay and molecular 
tumbling. Evidently these mechanisms can account for 
all relaxation experimentally observed in nl In HCA I 
perturbed angular correlation spectra. Since the coordi- 
nation complex initial metastable state is assumed to be 
independent of the Auger emission intensity, the chemi- 
cal transmutation model predicts equality of the K X ray 
coincidence and anti-coincidence spectra. 

Relaxation in the dose dependent damage model spec- 
tra results from ligand reorientation. The absence of 
molecular tumbling relaxation in this model is a matter of 
convenience. Additional high precision lllm Cd HCAB I 
spectra would be required to establish the correct balance 
of relaxation due to internal and overall motion. The the- 
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TABLE I: Chemical transmutation model quadrupole inter- 
action parameters for perturbed angular correlation spectra 
of HCA I. 



l Cd 



Metastable state 

Freq. mean (s _1 ) 

Freq. SD 

Asymmetry 

Lifetime (ns) 
Ground state 

Freq. mean (s" 1 ) 

Freq. SD 

Asymmetry 

Euler angles (<f>,d,ip) 
Tumbling lifetime(ns) 



0.0 

95 x 10 6 

18% 

0.68 

81 



0.08 



0.06 



0.04 



0.02 



133 x 10 b 

50% 

0.0 

66 

95 x 10 6 

18% 

0.68 

0,7r/2,0 
81 




20 40 60 80 100 120 

t (ns) 



TABLE II: Dose dependent damage model quadrupole inter- 
action parameters. 



Freq. mean (s x ) 
Freq. SD 
Asymmetry 
Rotation SD 
Jump lifetime (ns) 



electrons (R. W. Howell, private communication). The 
anti-coincidence yield must be corrected for the efficiency 
of detection for K X rays. The corrected anti-coincidence 
yield. 



Ulm Cd 


ln In (k+) 


ni In (k- 


95 x 10 6 


127 x 10 6 


133 x 10 6 


18% 


47.6% 


53% 


0.68 


0.468 


0.43 


25° 


37.7° 


40° 


26.5 


18.1 


16.5 



_ (l-") y t | "(1-6) 

1 — eoj 1 — euj 



k+, 



(3) 



where e is the K X ray detection efficiency, uj is the prob- 
ability per decay of K X ray emission, Yfc_ is the perfect 
anti-coincidence yield, and Y^ + the K X ray coincidence 
yield. At the estimated detection efficiency for cadmium 
K X rays of 23%, the corrected anti-coincidence Auger 
electron yield is 13.0. The ratio Y^/Y^ = 11/13 in- 
terpolates the K X ray coincidence parameters between 
the zero dose and K X ray anti-coincidence parameters 
in the dose dependent damage model. The ratio of the 
difference between the K X ray coincidence and zero dose 
parameters to the difference between the K X ray anti- 
coincidence and zero dose parameters is 11/13, see Ta- 
ble [n] The experimental 111 In HCA I K X ray coinci- 
dence and anti-coincidence spectra are identical. How- 
ever, since the difference in the theoretical K X ray co- 
incidence and anti-coincidence spectra is small enough 
relative to the experimental error, the dose dependent 
damage model may not be confidently rejected. 



FIG. 3: Dose dependent damage model for HCA I perturbed 
angular correlation spectra. The experimental data is identi- 
cal to that displayed in Fig. 2. The solid curve is an approx- 
imate fit of the dose dependent damage model to the data. 
The dotted curve is an approximate fit to the lllm Cd HCA I 
perturbed angular correlation spectrum with relaxation from 
ligand reorientation rather than molecular tumbling. The 
dashed curve is the dose dependent damage model for the K 
X ray coincidence HCA I perturbed angular correlation spec- 
trum derived by linear interpolation of the model parameters 
for the lllm Cd and K X ray anti-coincidence spectra. 

oretical K X ray coincidence and anti-coincidence spectra 
for the dose dependent damage model depend upon the 
relative Auger emission intensities. The electron cap- 
ture decay of 11 Hn in coincidence with a cadmium K 
X ray yields on average 11.0 Auger and Coster-Kronig 
electrons and in anti-coincidence yields on average 16.8 



Conclusions 

The indistinguishability of the experimental m In 
HCA I K X ray coincidence and anti-coincidence spec- 
tra suggests that during the first 10 -8 seconds following 
electron capture decay the hot atom chemistry is inde- 
pendent of the Auger and Coster-Kronig electron emis- 
sion intensity. As demonstrated by the dose dependent 
damage model, the experimental data do not conclusively 
rule out dependence of hot atom chemistry on emission 
intensity. The difference between the HCA I K X ray 
coincidence and anti-coincidence spectra could be signif- 
icantly increased by detecting X rays for the coinci- 
dence spectrum, increasing the detection efficiency of all 
K X rays for the anti-coincidence spectrum, and possibly 
by measuring at liquid nitrogen temperature. 

Special thanks are due to Kandula S. R. Sastry and 



6 



Franklyn G. Prendergast for their continued encourage- 
ment and support of this project. Lynda McDowell pu- 
rified HCA I and helped with indium binding assays. 
Heiner Winkler introduced me to the Monte Carlo proce- 
dure [2l| for computing a perturbation factor by averag- 
ing over the probability distribution of piecewise static 
interactions. Roger Howell provided electron yields. 
Alexander Halpern has called to my attention refer- 
ences 0, 0, IE El • This work was supported in part by 
National Institutes of Health grant 1 R03 RR02219-01. 



[1] 
[2] 



[3] 

M 
[5] 

[6] 
[7] 



Electronic address: haydock@mailaps.org Current ad- 
dress 1326 2 nd St. NW,' Rochester, MN 55901 
J. P. Adloff, Radiochim. Acta 25, 57 (1978). 
P. Boyer and A. Baudry, in Hot Atom Chemistry, edited 
by T. Matsuura (Elsevier, Amsterdam, 1984), pp. 315- 
347. 

H. Sano and P. Giitlich, in Hot Atom Chemistry, edited 
by T. Matsuura (Elsevier, Amsterdam, 1984), pp. 265- 
302. 

M. Alflen, C. Hennen, F. Tuczek, H. Spiering, P. Giitlich, 
and Z. Kajcsos, Hyperfme Interact. 47, 115 (1989). 
C. Haydock and K. S. R. Sastry, in Proceedings of the 
8th International Congress of Radiation Research (1987), 
vol. 1, p. 67. 

T. Kobayashi, K. Fukumura, T. Kitahara, and 
S. Shimizu, J. Phys. Colloq. pp. C28-29 (1979). 
C. M. Lederer and V. S. Shirley, eds., Table of Isotopes 
(John Wiley, New York, 1978), pp. 516-523, seventh ed., 



[8 
[9 

[10 

[11 

[12 
[13 



[14 

[is; 



[16 
[17 



[18 
[19 



[20 
[21 



A = 111, URL http: //ie . lbl .gov/toi/nucSearch. asp 
H. Kim and D. D. Dlott, J. Chem. Phys. 94, 8203 (1991). 
B. Diehn, A. Halpern, and G. Stocklin, J. Am. Chem. 
Soc. 98, 1077 (1976). 

R. G. Khalifah, D. J. Strader, S. H. Bryant, and S. M. 
Gibson, Biochemistry 16, 2241 (1977). 
J. B. Hunt, M.-J. Rhee, and C. B. Storm, Anal. Biochem. 
79, 614 (1977). 

Y. Pocker and J. T. Stone, Biochemistry 6, 668 (1967). 
H. Frauenfelder and R. M. Steffen, in Alpha-, Beta- 
and Gamma- Ray Spectroscopy, edited by K. Siegbahn 
(North- Holland, Amsterdam, 1965), vol. 2, chap. XIXA, 
pp. 997-1198. 

E. Matthais, W. Schneider, and R. M. Steffen, Ark. Fys. 
24, 97 (1963). 

H. Goldstein, Classical Mechanics (Addison- Wesley, 
Massachusetts, 1980), pp. 606-610, 2nd ed., appendix 
B. 

R. Bauer, Q. Rev. Biophys. 18, 1 (1985). 
B. T. Smith, J. M. Boyle, J. J. Dongarra, B. S. Gar- 
bow, Y. Ikebe, V. C. Klema, and C. B. Moler, Ma- 
trix Eigensystem Routines - EISPACK Guide (Springer- 
Verlag, Berlin, 1976), no. 6 in Lecture Notes in Computer 
Science, 2nd ed. 

F. James, Rep. Prog. Phys. 43, 1145 (1980). 

M. P. Allen and D. J. Tildesley, Computer Simulation of 
Liquids (Clarendon Press, Oxford, 1989), pp. 345-351, 
appendix G. 

R. Bauer, P. Limkilde, and J. T. Johansen, Biochemistry 
15, 334 (1976). 

E. Gerdau, H. Winkler, B. Giese, W. Gebert, and 
J. Braunsfurth, Hyperfme Interact. 1, 469 (1976). 



